arXiv:1507.01979v2 [cond-mat.str-el] 9 Nov 2015 


Influence of interstitial Mn on local structure and magnetism in Mnj^sSb 


Joshua A. Kurzman, Andrew J. Martinolich, and James R. NeilsoiQ 
Department of Chemistry, Colorado State University, Fort Collins, Colorado 80523-1872, United States 

(Dated: November 11, 2015) 

We report x-ray total scattering and pair distribution function (PDF) studies of the structural relaxation around 
interstitial manganese (Mn,) in ferromagnetic Mnj^gSb (0.03 < 5 < 0.23) alloys, guided by density functional 
theory (DFT). Refinements to the experimental PDF using a crystallographically constrained structural model 
indicate an expansion in the equatorial plane of the Mn/Sbs trigonal bipyramidal site, which introduces signif¬ 
icant positional disorder in addition to the nominally-random occupation of interstitial voids. Observation of 
a weak diffuse signal near the symmetry-forbidden ( 001 ) reflection position is indicative of correlated disor¬ 
der from the clustering of Mn,-. Density functional relaxation of supercells approximating the 5 = 0.08, 0.15, 
and 0.23 compositions provides improved models that accurately describe the short-range structural distortions 
captured in the PDFs. Such structural relaxation increases the DFT calculated moment on Mn,-, which aligns 
antiparallel to the primary Mn moments, but leads to insubstantial changes in the average Mn and Sb moments 
and moments of Mn and Sb proximal to interstitials, thus providing a more accurate description of the observed 
bulk magnetic properties. 


I. INTRODUCTION 


Manganese pnictides of the hexagonal NiAs structure type 
{e.g. MnAs, Mni_|_ 5 Sb, MnBi; Figure display a wealth 
of functional magnetic properties. The class of materials 
includes potential candidates for rare earth-free permanent 
magnets, DEI mag neto-optical data storage,l250 and magnetic 
refrigeration.EElManganese antimoni de is known to only exist 
as a manganese-rich phase, Mnj_|_gSb.l^li21 The wide composi¬ 
tional region shown on the equilibrium phase diagram extends 
from 0.02 < 5 < 0.23,^ with excess Mn occupying the in¬ 
terstitial site (Figure [^) in trigonal bipyramidal coordination 
by Sb. Increasing amounts of interstitial Mn (Mn,) lengthen 
the a and shorten the c lattice parameters, which accompanies 
reductions of the saturation magnetization, spin-reorienta tion 
temperature ( Tsr ), 22011 and the Curie temperature ( 7 c ).®3 
Magnetic ordering of the primary Mn sites occurs parallel to 
c at high temperature, but it changes to the ab plane below 
7sr.22 

Although the bulk magnetic properties are well known, the 
microscopic nature of the magnetism associated with Mn, 
has been equivocal. A straightforward explanation of the re¬ 
duced magnetization in more Mn-rich samples is that Mn,- is 
aligned antiparallel to the ferromagnetic MnMn atoms, that is, 
Mni_|_ 5 Sb is ferrimagnetic. (Here we use Kroger-Vink nota¬ 
tion to differentiate interstitial Mn, and nominal MnMn atoms.) 
Nonpolarized neutron diffraction studies are consistent with 
this expectation.22S2 However, polarized neutron scattering 
analyses have suggested that Mn, has no moment.^SSni The 
reduction of magnetization is then explained on the basis of 
greater orbital overlap between Mn atoms upon reduction of 
the c lattice constant, and to a local perturbation of the Mn 
moment when in proximity to Mn, .223El] 

The question of whether Mn, carries a moment was recently 
re-addressed by Taylor et al. via a combination of nonpolar¬ 
ized elastic and inelastic neutron scattering on an Mni i 3 Sb 
single crystal.22 A critical feature of the modeling strategies 
adopted in previous polarized neutron studies wa s the use of 
a highly aspherical magnetic form factor for Mn.lSEH! in the 
recent work of Taylor and co-workers it was shown that such 
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FIG. 1. (a) NiAs structure exemplified by the hard ferromagnetic 
phase MnBi [space group F 63 / mmc, with Mn on 2a at (0,0,0) and Bi 
on 2 c at (j, j, |)], and inelastic neutron scattering data collected at 
4 K from a single crystal of MnBi. (b) Structure of Mni_|_gSb high¬ 
lighting the fractionally occupied, trigonal bipyramidal coordination 
environment of interstitial Mn [Mn,-, on 2d site at ( 3 , |)]. Note 

that Mn; is also in trigonal prismatic coordination with respect to the 
primary (fully occupied) MnMn site. Inelastic neutron scattering data 
collected at 10 K on a single crystal of Mni i 3 Sb shows a pronounced 
diffuse component that is not observed for MnBi. Inelastic scattering 
figures adapted from Taylor et al., referencell4l 


a model provides an inferior description of the magnetic re¬ 
flections - particularly at low Q where magnetic scattering is 
strongest - relative to a model employing a spherical magnetic 
form factor and antiferromagnetic coupling between Mn, and 
MnMn -22 Density functional calculations enlisting small su¬ 
percells further supported Mni_,_gSb as a ferrimagnet. Addi¬ 
tionally, a pronounced and gapless diffuse magnetic compo¬ 
nent was observed across all temperature regimes in the in¬ 
elastic spectrum (Figure [^), centered at (001). Simil ar or 
related features have also been noted by other groups.23211 
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The location of the diffuse component is intriguing because 
the (001) reflection should be systematically absent from the 
structure factor; there is no diffuse scattering observed for sto¬ 
ichiometric (interstitial-free) MnBi (Figure [^). Two plausi¬ 
ble explanations were suggested; the diffuse scattering could 
arise from (/) correlated structural or magnetic disorder, or (ii) 
modification of the neighboring Mnjvin moments by Mn,P3 
An essential feature of the Mni_|_6Sb system that has been 
largely unaccounted for, both experimentally and computa¬ 
tionally, is the local structure around Mn,-. It has long been 
recognized that the equatorial Mn,—Sb contacts of the crys¬ 
tallographic Mn, site in the Mni_|_ 5 Sb unit cell are quite 
shortHowever, few studi es hav e attempted to model the 
presumed structural relaxationfSEl! Yamaguchi and Watanabe 
approximated that Mn, would be accommodated by distor¬ 
tions of only the Sb positions, the magnitude of the local dis¬ 
placements being derived from the concentration dependence 
of the lattice parameters!^ This assumption was used in their 
analysis of polarized neutron diffuse scattering data. Mod¬ 
els with no Mn, moment, but containing reduced moments on 
the first or second nearest-neighbor Mn atoms in proximity to 
Mn,, were compared against calculated scattering for a model 
with Mn,- aligned antiparallel. Although the best agreement 
was obtained for a reduction of the six nearest-neighbor Mn 
moments, an aspherical Mn magn etic form factor appears to 
have contributed to the result.E^lSl 

The utility of density functional modeling in describing lo¬ 
cal structural effects and defect physics is well demonstrated. 
A variety of approaches can be used to simulate defects (sub¬ 
stitutional, interstitial, antisite, vacancy) and random alloying, 
the validity of which is reflected in agreement with a vari¬ 
ety of experimental probes of local bonding arrangements, in¬ 
cluding nuclear magnetic resonance (NMR)|— extended x- 
ray absorption fine structure (EXAFS),l223Eni and pair distribu¬ 
tion function (PDF) data obtained from total scattering.l^DEI] 
When the concentration of defects is low, an effective strat¬ 
egy is to model a defect in a sufficiently large supercell such 
that the defect is approximately isolated (i.e., limiting the in¬ 
teraction between defects when periodic boundary conditions 
are applied).® At higher levels of substitution, the creation 
of models on the basis of chemically reasonable, ordered su¬ 
percells c an be an elegant strategy to assess plausible bond¬ 
ing motifs.®® Random alloys are inherently challenging sys¬ 
tems to model in the framework of density functional the¬ 
ory (DFT), since even large supercells will necessarily con¬ 
tain elements of periodicity that are absent in a truly random 
structure, but strategies such as stochastic mixing and the spe¬ 
cial quasirando m stru ctures approach often provide adequate 
approximations.®®ES!Mnj^gSb can be regarded as a type of 
random alloy, with a distribution of interstitials and vacancies 
over the crystallographic Mn, site. 

In the present contribution, we examine local structural re¬ 
laxation in Mnj^gSb using synchrotron x-ray total scattering 
and PDF analysis coupled with density functional modeling. 
The use of x rays rather than neutrons allows us to base our 
analysis on the local nuclear structure without interference 
from magnetic scattering. Particular attention is devoted to 
the Mn,- environment and its calculated magnetic moment, as 


well as the compositional and structural influence on the mo¬ 
ment of ferromagnetically aligned MnMn- Structural relax¬ 
ation has a pronounced effect on the Mn, environment and 
enhances electronic localization by a lengthening of equato¬ 
rial Sb contacts, but it does not appear to influence the mo¬ 
ments on MnMn or Sb. A weak diffuse signal is observed near 
the symmetry-forbidden (001) Bragg position, suggestive of 
correlated disorder associated with clustering of Mn,-. These 
results are discussed in light of the recent work of Taylor et 
al.® 


II. MATERIALS AND METHODS 
A. Preparation of Mn i+gSb 

The series of polycrystalline Mni_|_gSb (0.01 < 5 < 0.23) 
samples were prepared in alumina crucibles from the ele¬ 
ments (Mn flake, NOAH Technologies Corporation, 99.99 %; 
Sb shot, Alfa Aesar, 99 %). Antimony was recrystallized be¬ 
fore use, and surface oxide was removed from Mn by heating 
overnight in an evacuated fused silica tube at 980 °C. Appro¬ 
priate stoichiometric mixtures sealed in evacuated silica tubes 
were melted at 930 °C for 16 h and then quenched in water. 
The obtained ingots were finely powdered in an agate mor¬ 
tar and pressed into pellets, then annealed under vacuum at 
700 °C for 48 h and water-quenched. 


B. Characterization 

Magnetization measurements were performed with a vi¬ 
brating sample magnetometer on a Quantum Design, Inc. Dy- 
nacool PPMS. Room-temperature x-ray powder diffraction 
and x-ray total scattering measurements were conducted at the 
Advanced Photon Source of Argonne National Laboratory. 

High-resolution synchrotron x-ray powder diffraction pat¬ 
terns were collected at the 11-BM-B beamline using an x-ray 
energy of about 30keV (k « 0.414 A).®® Samples were di¬ 
luted with amorphous Si 02 to reduce the effects of absorption. 
Riet veld a nalyses were performed within the GSAS/EXPGUI 

suite.™ 

X-ray total scattering measurements were performed at the 
11-ID-B beamline. Data were collected with an amorphous 
silicon area detectoi® using two x-ray energies, ^ 58keV 
(A. = 0.2114 A) and ^ 86keV (A, = 0.1430 A), at sample-to- 
detector distances of approximately 17 cm and 19 cm, respec¬ 
tively. Calibrations were performed by measurement of a 
Ce02 standard at each condition. The electron density pair 
distribution function, G(r), was obtained from background- 
subtracted scattering data by the ad-hoc approach applied in 
PDFgetX3.®The reduced scattering structure function, S{Q), 
was transformed to G{r) using a maximum momentum trans¬ 
fer of Q max — 24 A foi data collected at 58keV, and ^max — 
28 A^^ for data collected at 86keV. A powdered nickel stan¬ 
dard was used to determine the resolution truncation parame¬ 
ters Qiamp Abroad used in PDFgui refinements.® Reverse 
Monte Carlo (RMC) simulations were conducted using the 








3 


RMCprofile softwarej®! with models constrained to fit S{Q) 
and G{r) to capture both long- and short-range order. A “clos¬ 
est approach” constraint of 2.2 A was applied to all atomic 
species to prevent the fitting of termination ripples in the low 
r region of the PDF. Simulated diffraction patterns were com¬ 
puted using the DIFFaX software® input files are provided 
in the supplemental material. 


C. Computational Details 

Density functional theory (DFT) calculations were per¬ 
formed using the Vienna ab-initio simulation package 
(VASP),® with interactions between the cores (Mn:[Ar], 
Sb:[Kr]4d^*’) and valence electrons described using the pro¬ 
jector augmented wave (PAW) method.® Calculations were 
performed within the generalized gradient approximation 
(GGA.) using the functional of Perdew, Burke, and Ernzerhof 
(PBE)^ to account for the effects of exchange and correla¬ 
tion. An energy cutoff of 350 eV and a F-centered 4x4x4 
k-point sampling were employed for all calculations. Since 
PDF analysis is a measurement of the ensemble average of 
all pairwise correlations in a sample, a series of 10 dis¬ 
tinct initial supercell configurations were constructed for each 
of three compositions, Mn 48 (Mn;) 4 Sb 48 , Mn 48 (Mn,) 7 Sb 48 , 
and Mn 48 (Mn,)iiSb 48 (where Mn, denotes interstitial man¬ 
ganese), which respectively provide compositions close to 
Mni.osSb, Mni.isSb, and Mni 23 Sb. Parent Mn48Sb48 super¬ 
cells were built from orthonormal C-centered orthorhombic 
supercells of the conventional crystallographic lattice, with di¬ 
mensions reflecting the lattice constants obtained by Rietveld 
refinement. Interstitial Mn positions were selected stochasti¬ 
cally with the aid of RMCprofild^ by filling all of the inter¬ 
stitial positions with Mn or vacancy (dummy) atoms and dis¬ 
tributing Mn,- by performing short simulations in which Mn, 
atoms were “swapped” with vacancy atoms; simulations last¬ 
ing 6 seconds typically produced 50,000 atom swap moves, 
leading to a robust ensemble of starting structural configura¬ 
tions. Relaxations of the supercells were performed at fixed 
cell dimensions, allowing only the atomic positions to change, 
and they were deemed to have converged when the forces on 
all ions were less than 0.01 eVA^*. All calculations were 
performed with spin polarization, using starting spin config¬ 
urations in which Mn, was aligned antiparallel to Mn. For 
select configurations, test calculations were performed with 
an on-site Coulomb repulsion term on Mn (GGA+U), or with 
spin-orbit coupling; neither the inclusion of a Hubbard U nor 
the effects of the spin-orbit interaction produced any qualita¬ 
tive changes to the results discussed below. 
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FIG. 2. (a) Axial (c/a) ratio of the Mnj^gSb lattice constants as a 
function of the nominal interstitial Mn content, 5, showing Vegard 
law behavior over most of the compositional range and excellent 
agreement between Rietveld and PDF refinements and literature 
values.l3 Error bars are smaller than the symbols, (b) Room tem¬ 
perature saturation magnetization as a function of 8, and bulk mag¬ 
netization at 0 K obtained from DFT calculations of unrelaxed and 
relaxed supercells of Mn 48 (Mn,)j:Sb 48 (x = 4,7,11); DFT moments 
represent the average of 10 distinct configurations at each composi¬ 
tion, with error bars representing the standard deviation of the av¬ 
erage. (c) Room temperature magnetization hysteresis curves of the 
Mnj^gSb series of samples. Note that the magnetization curves for 
8 = 0.01 and 8 = 0.03 nearly overlay, which is consistent with the 
samples’ c/a ratios being nearly the same. 


III. RESULTS AND DISCUSSION 

A. Crystallographic analyses and bnlk magnetic properties 

The series of Mnj^gSb samples (0.01 < 5 < 0.23) span 
approximately the entire compositional range and display 
Vegard law behavior of the lattice parameters between 0.03 < 


5 < 0.23, as shown by the linear trend in c/a ratio as a function 
of 5, Figure [^. Small deviations from linearity are noted for 
the nominally 5 = 0.01 and 5 = 0.23 compositions, indicating 
proximity to the phase boundary extrema. The 5 = 0.01 sam¬ 
ple contains a small amount of antimony metal (< 1 wt. %), 
and all of the samples contain small impurities of MnO, esti¬ 
mated by quantitative phase analysis to be < 1 wt. %. There 
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is no evidence of an Mn 2 Sb impurity in the Mni 23 Sb sample, 
suggesting the solubility limit had not been reached. Axial 
ratios obtained from Rietveld refinements of high-resolution 
synchrotron x-ray diffraction data (11-BM) agree very closely 
with the ratios obtained from real-space refinements of PDF 
data. A representative Rietveld refinement is presented in Fig¬ 
ure]^ for the 5 = 0.05 sample, and a complete list of refined 
values for both Rietveld and PDF refinements is given in Ta¬ 
ble SM-1. The Mn, occupancies estimated by Rietveld refine¬ 
ment agree well with the nominal stoichiometries across the 
0.03 < 5 < 0.23 range, with estimated errors in the site occu¬ 
pancy of 0.01. Fractionally occupied Mn; has a site multiplic¬ 
ity of 2 (Wyckoff position 2d) in the crystallographic model, 
but this does not imply that interstitial atoms necessarily oc¬ 
cur in parrs in a given unit cell of the “real” material. Consider 
that for 5 = 0.05 there is only one Mn, for every 10 unit cells 
[i.e. Mn 2 o(Mn,)iSb 2 o]: although the crystallographic cell in¬ 
dicates that Mn, populates both sites of the 2d position on av¬ 
erage, it is quite conceivable that interstitial atoms are present 
locally in just one of the two interstitial voids. 
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FIG. 3. Representative Rietveld refinement of Mn; osSb against high 
resolution synchrotron x-ray diffraction data (11-BM, Advanced 
Photon Source). To aid in visual inspection of the fit, the low and 
high Q regions of the diffraction data are omitted. Reflection posi¬ 
tions for Mni osSb and the MnO impurity are shown as solid vertical 
lines in the two middle panels. 

Room-temperature saturation magnetizations (Figure |^), 
taken from the room-temperature hysteresis curves shown in 
Figure |^, decrease linearly as the quantity of interstitial Mn 
increases. Magnetization curves of the 5 = 0.01 and 5 = 0.03 
samples are nearly overlaid, which is consistent with their 
comparable c/a ratios and 5 values obtained by Rietveld re¬ 
finement. 

In general, PDFs of the Mnj^gSb series are well de¬ 
scribed by the conventional crystallographic cell (constrained 
by space-group symmetry), particularly at low concentrations 
of Mn, where there is limited sensitivity to its pair cotTela- 
tions. Real-space refinements (PDFgui) over the 2-30 A range 
for Mni.oaSb, Mn; i 2 Sb, and Mn; 23 Sb are shown in Figure 


and they are representative of fit quality across the com¬ 
plete series of samples (Figure SM-1). These fits worsen as a 
function of increasing 5 in the low r region of the PDFs, but 
beyond about 5 A the crystallographic structures provide very 
satisfactory fits to the data. In fact, for samples containing up 
to 10% Mn,-, better statistical agreement is consistently ob¬ 
tained in refinements over the 6 to 30 A range with a Mn; poSb 
model (i.e., without Mn,) than for Mni_|_gSb models contain¬ 
ing partially occupied Mn,. 



FIG. 4. PDF refinements (solid lines) against data (open circles) 
for samples of nominal compositions Mn; gsSb, Mnj i 2 Sb, and 
Mni 23 Sb, employing the crystallographic structural model in which 
interstitial Mn (Mn,) is located on the 2d site and partially occu¬ 
pied. (a) Full-range fits, from 2 to 30 A. The vertical line at r = 6 A 
demarcates different x-axis scales at shorter and longer r. (b) Pan¬ 
els showing the first-nearest-neighbor region of the three full-range 
fits; select PDF partial contributions for MnMn-Sb (dotted line) and 
Mn,-Sb (solid line) correlations are vertically offset for clarity. As 
the amount of Mn, increases, the crystallographic structural model’s 
deficiency becomes more apparent: there are shorter pair correlations 
present in the model than observed in the data, and these arise from 
equatorial Mn;-Sb distances. 
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The corresponding first nearest-neighbor peaks, compris¬ 
ing contributions from MnMn-Sb, Mn,-Sb and Mn,—MnMn, 
are highlighted in Figure [^; the PDF partials for Mumh- 
Sb and Mn,—Sb pair correlations are shown below in dotted 
and solid lines, respectively. For clarity, MnMn-Mn, partials 
are not shown as these contribute very weakly at the same 
distances as MnMn-Sb. The crystallographic model predicts 
contributions at approximately 2.4 A, although none are ob¬ 
served experimentally; this intensity is derived from equato¬ 
rial Mn,-Sb correlations of the Mn, trigonal bipyramid. The 
absence of observed intensity at this distance in G(r) is ev¬ 
idence that structural relaxation occurs either by or around 
interstitial atoms, which is also reflected in the poorer de¬ 
scription of peak intensity and breadth for the more Mn-rich 
samples. This result is quantitatively consistent for the PDFs 
obtained from total scattering with two different x-ray wave¬ 
lengths. This feature is also robust considering the Nyquist 
sampling frequency dictated by Qmax (Figure SM-2).I^ Com¬ 
paring the crystallographic structural model against the ob¬ 
served data indicates that the equatorial Mn,-Sb contacts are 
longer in the “real” material than reflected by the crystallo¬ 
graphic cell. Qualitatively, a bond-lengthening around Mn, 
would be expected to increase electron localization and its 
magnetic moment.El The ability to directly inspect PDF data 
and infer plausible structure-property relations is an indis¬ 
putable asset of the technique. However, a more rigorous 
modeling approach is required to support and reconcile the 
impact of structural relaxation on the properties of Mni_|_ 6 Sb. 


B. DFT: unifying local structure and magnetic properties 

A series of density functional calculations were performed 
to assess relaxation by and around interstitial atoms in 
Mni_|_ 5 Sb. To account for the disorder inherent to this sys¬ 
tem, a collection of 10 distinct supercell configurations was 
created for each of three compositions selected from the com¬ 
positional range. All configurations are based on 96 atom 
MnSb supercells, with Mn 48 (Mn,) 4 Sb 48 , Mn48(Mn,)7Sb48, 
and Mn 48 (Mn,)iiSb 48 formulas chosen to provide composi¬ 
tions close to Mni ogSb, Mni.isSb, and Mni 2386 . The po¬ 
sitions of Mn, in the supercell configurations were selected 
stochastically. Representative structural depictions of one 
configuration for each composition are shown in Figure]^ To 
facilitate comparison with the PDF data, supercells for a given 
composition were constrained to dimensions commensurate 
with the experimentally determined unit-cell parameters. 

All of the PDFs for different configurations of a given 
composition are very similar (Figure SM-3), and differences 
are far less significant than termination artifacts introduced 
in the Fourier transform of S{Q). A comparison of fits at 
short r to the PDF of Mni 2386 for the crystallographically- 
constrained structure (unit-cell model) versus a representa¬ 
tive Mn 48 (Mn,)iiSb 48 supercell is shown in Figure]^. DFT 
relaxation provides a clear improvement in the description 
of G(r), reflected in the marked reduction of R„', for fits 
over the range 2 to 6 A the average R-factor for all 10 con¬ 
figurations is 13(1)%, whereas R„ = 24% for the crystal- 



Mn,g(Mn,,)„Sb,3 


FIG. 5. Structural depictions of representative relaxed supercell con¬ 
figurations for Mn48(Mn,)4Sb48 (top), Mn 48 (Mn,) 7 Sb 48 (middle), 
and Mn 48 (Mn,)iiSb 48 (bottom). MnMn is shown as light (pink) 
spheres, Sb as dark (gray) spheres, and Mn, (blue) in polyhedral ren¬ 
dering. An outline of the conventional unit cell is shown for refer¬ 
ence. 


lographic structural model. Statistical improvement with 
the Mn 48 (Mn,) 7 Sb 48 DFT relaxed structures is also found 
for fits against the Mni 1586 PDF data (not shown), with 
R„ = 13(1)% versus R„ = 16%. On the other hand, re¬ 
laxed Mn 48 (Mn,) 4 Sb 48 configurations are marginally worse 
than the crystallographic structure, statistically, in fits against 
the Mni,o 88 b data (R„ = 17.1(4) % vs. R„ = 15.1 %), but this 
most likely reflects limited sensitivity to Mn, correlations for 
this composition. Indeed, refinement of a Mni 008 b unit-cell 
model against the Mni 0386 data provides a better statistical 
fit (R„ = 13.5 %) than such a model with partially occupied 
Mn,. 

Select atom-pair partials for the crystallographic Mni 238 b 
structure and relaxed Mn 48 (Mn,) 118648 supercell are shown in 
Figure |^. The only pronounced differences occur for Mn,- 
Sb and Sb-Sb pair correlations, with MnMn- 8 b and MnMn- 
MnMn partials being notably very similar for the relaxed and 
unit cell models. As anticipated, equatorial Mn,-Sb correla- 
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FIG. 6. (a) Comparison of PDF refinements for a 

crystallographically-constrained (unit cell) structural model and 
representative DFT-relaxed configuration against data collected for 
the sample of nominal composition Mni 23 Sb. Refinement R-factors 
shown from fits in the range of 2 to 6 A; Rw for the DFT relaxed 
configuration represents the average for all 10 configurations, given 
with the standard deviation of the average, (b) Select atom-pair 
partials extracted from the DFT-relaxed and unit cell models, along 
with PDF partials generated from a reverse Monte Carlo simulation 
using a 9,990 atom supercell {R„ = 8.6 %). 


tions are clearly shifted to longer r in the DFT-relaxed model. 
Differences in the Sb-Sb correlations of the models indicates 
that the antimony sublattice experiences the largest disruption 
upon incorporation of Mn, . PDF partials for Mn,-Sb and Sb- 
Sb correlations obtained from reverse Monte Carlo simulation 
(bottom of Figure]^) agree well with those obtained by DFT, 
lending additional support to our approach. 

Histograms displaying the amplitudes of atom displace¬ 
ments away from their crystallographic (average) structural 
positions upon DFT relaxation are shown in Figure ar¬ 
ranged by site, and combined for all 10 supercell configura¬ 
tions of each composition. The Mnjvin and Mn,- sites show uni- 
modal distributions extending to roughly 0.1, 0.15, and 0.2 A, 
respectively, from lowest to highest Mn, content. The approx¬ 
imately bimodal Sb histograms show displacements about 
twice as large for a given composition, with a third region of 
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FIG. 7. Histograms showing the amplitude of atom displacements 
away from the crystallographic structural positions upon DFT relax¬ 
ation, compiled for all 10 configurations of each DFT supercell com¬ 
position. Top row: Mn 4 g(Mn/) 4 Sb 48 ; middle row: Mn 4 g(Mn,) 7 Sb 48 ; 
bottom row: Mn 48 (Mn,) 11 Sb 48 . 


lower frequency but larger amplitude displacements observed 
for the Mn 48 (Mn,) 7 Sb 48 and Mn 48 (Mn,)iiSb 48 configurations 
associated with Sb that have more than one nearest-neighbor 
Mn,. All of the larger amplitude Sb displacements (> 0.1 A) 
correspond to equatorial Mn,-Sb distances, whereas the dis¬ 
placement of axial Sb atoms is typically < 0.1 A. Despite the 
large degree of relaxation displayed by Sb atoms, the local 
distortions necessary to accommodate Mn, induce very mini¬ 
mal changes in the Mnjvin-Sb and MnMn-MnMn bond distance 
distributions (Figure]^), which are preserved because of co¬ 
operative distortions by MnMn- This suggests that relaxation 
should have a much more pronounced effect on the electronic 
environment of Mn,- than on either MnMn or Sb. 

The average calculated moments of the MnMn, Sb, and Mn, 
sites in the DFT relaxed supercells, along with the correspond¬ 
ing moments of unrelaxed supercells, are shown in Figure]^. 
All 10 configurations are shown for each composition, but 
there is no relation between configurations of the same num¬ 
ber for different compositions. It is immediately apparent that 
the moments on MnMn and Sb are insensitive and effectively 
invariant with respect to the different compositions, the differ¬ 
ent configurations, and also to structural relaxation. For the 
MnMn and Sb sites, the standard deviations of the average mo¬ 
ments are much smaller than the symbols. The average mo¬ 
ment on Mn,, however, is consistently of smaller magnitude 
in unrelaxed configurations than the average Mn, moment in 
relaxed supercells, but it is generally opposite in sign to that 
of MnMn- This is reflected in a smaller net magnetization for 
relaxed structures (Figure [^). For many of the supercells, 
including both relaxed and unrelaxed structures, there is con¬ 
siderable variability in the individual moments of Mn, sites, 
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Mn48(Mn.)4Sb48 Mn43(Mn,)^Sb4g Mn4g(Mn,)^^Sb43 



FIG. 8. (a) Average moment, by site, for each of the relaxed (filled symbols) and unrelaxed (open symbols) configurations among the three 
compositions examined by DFT. Error bars, representing one standard deviation of the average moments, are much smaller than the symbols for 
MnMn and Sb. Note that the average moment on MnMn and Sb are quite insensitive to structural relaxation, whereas the magnitude of the mo¬ 
ment on Mn,' becomes significantly larger when accounting for relaxation of and around the interstitials, (b) Energy of the relaxed (unrelaxed) 
configurations relative to the most stable relaxed (unrelaxed) configuration for a given composition. The energy stabilization of structural 
relaxation is about 20meV atom^* for Mn 48 (Mn,) 4 Sb 4 g, 30meV atom^* for Mn 48 (Mn/) 7 Sb 4 g, and 40meV atom^* for Mn 48 (Mn,)iiSb 48 . 


indicated by the error bars which show the standard devia¬ 
tion of the averages. The individual moments for each Mn, 
site in each configuration are shown in Figure SM-4 of the 
supplemental material. The configuration coordinates are also 
supplied therein. 

Types of variability in the Mn, moments found in individual 
configurations include the following: (a) sites in relaxed con¬ 
figurations that remain antiparallel but with decreased magni¬ 
tude of the moment, (b) sites in relaxed configurations with 
appreciable magnitude that display parallel alignment with 
respect to MnMn, (c) sites in unrelaxed configurations that 
are antiparallel and of significantly greater magnitude than 
the average moment, and (d) sites in unrelaxed configura¬ 
tions that have sizable moments and parallel alignment with 
MnMn. What is surprising about this result is that variabil¬ 
ity in the Mn,- moments does not necessarily appear to be as¬ 
sociated with an energetic penalty. The energies of the re¬ 
laxed configurations relative to the lowest energy (most sta¬ 
ble) relaxed configuration for a given composition are shown 
in Figure]^; the same comparison is also shown for unrelaxed 
configurations (relative to the most stable unrelaxed configu¬ 
ration). For example, configurations #9 are (coincidentally) 
the lowest energy structures for both Mn 48 (Mn,) 4 Sb 48 and 
Mn 48 (Mn,) 7 Sb 48 , including both the unrelaxed and relaxed 
supercells. Configuration 9 of Mn 48 (Mn,) 4 Sb 48 is seemingly 
well-behaved, with large antiparallel moments in the relaxed 
structure, and small antiparallel moments in the unrelaxed 


structure. In configuration 9 forMn 48 (Mn,) 7 Sb 48 , on the other 
hand, one of the Mn, sites in the relaxed structure has a sig¬ 
nificantly smaller (still antiparallel) moment, and two of the 
sites in the unrelaxed structure have notably larger moments 
than the average. 

Similarly, there are other configurations for each composi¬ 
tion that contain substantial variability in the individual Mn, 
moments but are lower energy structures than configurations 
with far less variability. In the case of the relaxed supercells, 
decreased moments are associated with shorter Mn,-MnMn 
and Mn,-Sb distances (Figure SM-5). This is consistent with 
the general observation that the Mn,- moments are smaller in 
unrelaxed supercells. It does not explain why some of the in¬ 
dividual moments are large in the unrelaxed structures; how¬ 
ever, it is conceivable that the observed variations result from 
the imposed periodicities. Nevertheless, the general trends in 
the MnMn, Sb, and Mn, moments should not be affected sig¬ 
nificantly. 


C. Diffuse scattering and correlated disorder of Mn, 

We now return to the question of the diffuse magnetic scat¬ 
tering observed at the nominally systematically absent (001) 
position in the recent report by Taylor and co-workers (refer¬ 
ence M- The signal was observed at the lowest energies ac¬ 
cessible in the inelastic experiment (~ 4 meV), but the pres- 
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ence or absence of an associated elastic diffuse signal could 
not be determinedfSI As noted by the authors, the observation 
of scattering intensity at this location - regardless of whether 
it is elastic or inelastic - implies symmetry breaking between 
the upper and lower portions of the unit cell along the c axis. 
Such breaking could be induced if some Mumh are antiferro- 
magnetically aligned along c (which would double the mag¬ 
netic unit cell), but both our work and theirs indicates that 
this is highly improbable. An alternative situation that would 
break the symmetry is the local population of Mn, in only one 
of the two interstitial voids of a unit cell, which removes the 
/ = 2n reflection condition. This latter possibility, where in¬ 
terstitial atoms do not necessarily occur in pairs in a given 
unit cell, was introduced earlier in the discussion of Rietveld 
refinements (Section III A). 


T-'- 1 -'- 1 -'- 1 -'-r 



FIG. 9. (a) Diffuse x-ray scattering features observed in the vicin¬ 
ity of the ( 001 ) reflection in Mnj^sSb samples (8 = 0 . 12 , 0.18, 
and 0.23); container scattering has been subtracted from the data, 
which was collected at ~ 86 keV at a sample-to-detector distance of 
~ 19 cm. (b) Simulated diffraction patterns of Mni 2386 structures 
built from an infinite number of infinitely-wide discrete layers. When 
layers with a singly occupied Mn, site per unit cell are preferentially 
stacked adjacent to layers with Mn, in the same site, the correlated 
disorder gives rise to a diffuse feature that is consistent with the ex¬ 
perimentally observed scattering in regards to both position (in Q) 
and relative intensity. Cartoon models of the different stmctures giv¬ 
ing rise to the simulated scattering are shown as four layer fragments; 
MnMn (pink), Sb (gray), Mn,- (blue). 


Close examination of our x-ray scattering data reveals a 
weak diffuse component on top of the systematically absent 
(001) reflection that is suggestive of correlated structural dis¬ 
order. In Figure we present the low-Q region of x-ray 
diffraction patterns for three Mni+gSb samples (5 = 0.12, 
0.18, and 0.23), from which the container scattering has been 
subtracted. Scattering intensities have been normalized to the 
most intense peak in the x-ray diffraction pattern of Mni^^gSb 
{i.e. the (101) reflection]. The intensity of the diffuse signal 
increases with 5, and the maxima of the diffuse features are 
all observed at lower Q (larger c/-spacing) than the center of 
the (001) reflection on the basis of the c lattice constant. A 
shift of diffuse scattering maxima away from the associated 
Bragg positions is commonly observed in situations of local 
structural dilation or contraction.El 

A series of simulations are presented in Figure]^ that offer 
support to the hypothesis that the diffuse signal at low Q arises 
from correlated structural disorder associated with “unpaired” 
Mn,, that is, occupation of just one of the interstitial sites in 
a unit cell. Diffraction patterns were generated with the DIF- 
FaX softward^ - which explicitly computes the incoherent 
intensity contribution to a given hkl reflection - for structures 
comprising random and correlated stacking variations of dis¬ 
crete layers. All simulations were conducted on the basis of 
the Mni 23Sb stoichiometry and cell dimensions. 

As a control DIFFaX simulation, a single-layer Mni_|_5Sb 
model containing partially occupied Mn, in both of the in¬ 
terstitial sites was considered. This model is designated “all 
Mn, paired” in Figure]^, and it results in no noticeable scat¬ 
tering intensity at the (001) position. Models containing un¬ 
paired Mn,- were constructed with three distinct layer types: 
(1) a layer without any interstitial atoms, (2) a layer with a 
fully occupied interstitial atom in the lower interstitial void, 
and (3) a layer with a fully occupied interstitial atom in the 
upper interstitial void. Two types of disorder were modeled 
for structures incorporating layers with singly occupied inter¬ 
stitial voids, namely random disorder and correlated disorder, 
as labeled in Figure For the model containing layers with 
unpaired Mn,- stacked in a completely random assortment, no 
resolved diffuse intensity is generated at the (001) reflection. 
The rise in intensity at low Q is the result of the diffuse scatter¬ 
ing from completely incoherent contributions, which emanate 
from (000). On the other hand, for an unpaired Mn, model in 
which there is preferential stacking of the layers - correlated 
disorder - a well resolved diffuse signal near the (001) is gen¬ 
erated. More specifically, the preferential stacking modeled 
here corresponds to an increased probability that a layer con¬ 
taining a single Mn, in one of the interstitial voids will neigh¬ 
bor layers that also contain Mn, in the same site. Although this 
simulation represents a rather coarse example of one type of 
correlated disorder that could be present in Mnj^gSb, it offers 
strong support for the presence of some degree of clustering 
associated with the interstitial atoms. As noted for the experi¬ 
mentally observed diffuse signal, the maximum of the diffuse 
component in the simulated pattern is also shifted to lower Q 
than the center of the (001) reflection. A single-crystal x-ray 
diffuse scattering study would be better suited to elucidating 
the origin of this observation. 
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While these observations and simulations soundly establish 
the existence of short-range structural order, the specific na¬ 
ture of clustering in Mni_|_5Sb remains an open question. Con¬ 
sideration of the many plausible local arrangements is beyond 
the scope of the present work; it is easily conceivable that a 
variety of short range-ordered motifs coexist. It is important 
to stress that although the presence of some unpaired Mn,- is 
a requisite for generating nuclear scattering intensity near the 
(001) reflection, this does not imply the absence of local re¬ 
gions with pairs of interstitials. Indeed, this type of diversity 
is represented in the DFT supercells. A comparison of the rel¬ 
ative supercell stabilities does not suggest an obvious penalty 
for interstitial proximity (pairing), but because of their rela¬ 
tively small sizes and construction by stochastic population of 
the interstitial sites it would be conjecture to suggest an ener¬ 
getic preference for any particular type of clustering. 

In summary, the observation of a diffuse signal consistent 
with the symmetry-forbidden (001) reflection suggests that 
the diffuse magnetic component observed by Taylor et al. by 
inelastic neutron scattering is rooted in short-range structural 
order. This explanation is further supported by our de nsity 
functional studies, which challenge previous notional^HIMlI 
that Mn, induces a significant reduction in neighboring MnMn 
moments. 


IV. CONCLUSIONS 

Analysis of x-ray total scattering data collected on 
polycrystalline powders of Mni_|_5Sb (0.03 < 5 < 0.23) 


reveals signatures of structural relaxation associated with 
the accommodation of interstitial atoms. In particular, the 
trigonal bipyramidal site of Mn, is better described by longer 
equatorial Mn,-Sb distances than reflected by the crystal¬ 
lographic unit cell. Density functional relaxation of large 
supercells for compositions close to 5 = 0.08, 0.15, and 0.23, 
provides excellent structural descriptions of the short-range 
disorder induced by Mn,, evidenced by improved fits to the 
experimental PDFs. Relaxation has a pronounced effect 
on the calculated moment of Mn,, which aligns antiparallel 
to MnMn. but our results suggest that Mn, has effectively 
no influence on the moments of MnMn in proximity to 
interstitials. The observation of a diffuse signal near the 
(001) reflection position, which is systematically absent when 
Mn, is distributed completely randomly, indicates that some 
degree of correlated disorder exists in Mni_|_gSb. 
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